Entanglement and Spontaneous Symmetry Breaking in Quantum Spin Models

It is shown that spontaneous symmetry breaking does not modify the ground-state entanglement of two spins, as defined by the concurrence, in the XXZ- and the transverse field Ising-chain. Correlation function inequalities, valid in any dimensions for these models, are presented outlining the regimes where entanglement is unaffected by spontaneous symmetry breaking

Reduction of the sign problem using the meron-cluster approach

The sign problem in quantum Monte Carlo calculations is analyzed using the meron-cluster solution. The concept of merons can be used to solve the sign problem for a limited class of models. Here we show that the method can be used to \textit{reduce} the sign problem in a wider class of models. We investigate how the meron solution evolves between a point in parameter space where it eliminates the sign problem and a point where it does not affect the sign problem at all. In this intermediate regime the merons can be used to reduce the sign problem. The average sign still decreases exponentially with system size and inverse temperature but with a different prefactor. The sign exhibits the slowest decrease in the vicinity of points where the meron-cluster solution eliminates the sign problem. We have used stochastic series expansion quantum Monte Carlo combined with the concept of directed loops.Comment: 8 pages, 9 figure

Field-induced XY behavior in the S=1/2 antiferromagnet on the square lattice

Making use of the quantum Monte Carlo method based on the worm algorithm, we study the thermodynamic behavior of the S=1/2 isotropic Heisenberg antiferromagnet on the square lattice in a uniform magnetic field varying from very small values up to the saturation value. The field is found to induce a Berezinskii-Kosterlitz-Thouless transition at a finite temperature, above which a genuine XY behavior in an extended temperature range is observed. The phase diagram of the system is drawn, and the thermodynamic behavior of the specific heat and of the uniform and staggered magnetization is discussed in sight of an experimental investigation of the field-induced XY behavior.Comment: 4 pages, 4 figure

Directed geometrical worm algorithm applied to the quantum rotor model

We discuss the implementation of a directed geometrical worm algorithm for the study of quantum link-current models. In this algorithm Monte Carlo updates are made through the biased reptation of a worm through the lattice. A directed algorithm is an algorithm where, during the construction of the worm, the probability for erasing the immediately preceding part of the worm, when adding a new part,is minimal. We introduce a simple numerical procedure for minimizing this probability. The procedure only depends on appropriately defined local probabilities and should be generally applicable. Furthermore we show how correlation functions, C(r,tau) can be straightforwardly obtained from the probability of a worm to reach a site (r,tau) away from its starting point independent of whether or not a directed version of the algorithm is used. Detailed analytical proofs of the validity of the Monte Carlo algorithms are presented for both the directed and un-directed geometrical worm algorithms. Results for auto-correlation times and Green functions are presented for the quantum rotor model.Comment: 11 pages, 9 figures, v2 : Additional results and data calculated at an incorrect chemical potential replaced. Conclusions unchange

Aspect-ratio dependence of the spin stiffness of a two-dimensional XY model

We calculate the superfluid stiffness of 2D lattice hard-core bosons at half-filling (equivalent to the S=1/2 XY-model) using the squared winding number quantum Monte Carlo estimator. For L_x x L_y lattices with aspect ratio L_x/L_y=R, and L_x,L_y -> infinity, we confirm the recent prediction [N. Prokof'ev and B.V. Svistunov, Phys. Rev. B 61, 11282 (1999)] that the finite-temperature stiffness parameters \rho^W_x and \rho^W_y determined from the winding number differ from each other and from the true superfluid density \rho_s. Formally, \rho^W_y -> \rho_s in the limit in which L_x -> infinity first and then L_y -> infinity. In practice we find that \rho^W_y converges exponentially to \rho_s for R>1. We also confirm that for 3D systems, \rho^W_x = \rho^W_y = \rho^W_z = \rho_s for any R. In addition, we determine the Kosterlitz-Thouless transition temperature to be T_KT/J=0.34303(8) for the 2D model.Comment: 7 pages, 8 figures, 1 table. Minor changes to published versio

Universal scaling at field-induced magnetic phase transitions

We study field-induced magnetic order in cubic lattices of dimers with antiferromagnetic Heisenberg interactions. The thermal critical exponents at the quantum phase transition from a spin liquid to a magnetically ordered phase are determined from Stochastic Series Expansion Quantum Monte Carlo simulations. These exponents are independent of the interdimer coupling ratios, and converge to the value obtained by considering the transition as a Bose-Einstein condensation of magnons, alpha_(BEC) = 1.5. The scaling results are of direct relevance to the spin-dimer systems TlCuCl_3 and KCuCl_3, and explain the broad range of exponents reported for field-induced ordering transitions.Comment: 4 pages, 4 eps-figure

Universal SSE algorithm for Heisenberg model and Bose Hubbard model with interaction

We propose universal SSE method for simulation of Heisenberg model with arbitrary spin and Bose Hubbard model with interaction. We report on the first calculations of soft-core bosons with interaction by the SSE method. Moreover we develop a simple procedure for increase efficiency of the algorithm. From calculation of integrated autocorrelation times we conclude that the method is efficient for both models and essentially eliminates the critical slowing down problem.Comment: 6 pages, 5 figure

Directed Loop Updates for Quantum Lattice Models

This article outlines how the quantum Monte Carlo directed loop update recently introduced can be applied to a wide class of quantum lattice models. Several models are considered: Spin-S XXZ models with longitudinal and transverse magnetic fields, boson models with two-body interactions, and 1D spinful fermion models. Expressions are given for the parameter regimes were very efficient "no-bounce" quantum Monte Carlo algorithms can be found.Comment: 18 pages, 19 figure

Two-Dimensional Quantum XY Model with Ring Exchange and External Field

We present the zero-temperature phase diagram of a square lattice quantum spin 1/2 XY model with four-site ring exchange in a uniform external magnetic field. Using quantum Monte Carlo techniques, we identify various quantum phase transitions between the XY-order, striped or valence bond solid, staggered Neel antiferromagnet and fully polarized ground states of the model. We find no evidence for a quantum spin liquid phase.Comment: 4 pages, 4 figure

Transverse Ising Model: Markovian evolution of classical and quantum correlations under decoherence

The transverse Ising Model (TIM) in one dimension is the simplest model which exhibits a quantum phase transition (QPT). Quantities related to quantum information theoretic measures like entanglement, quantum discord (QD) and fidelity are known to provide signatures of QPTs. The issue is less well explored when the quantum system is subjected to decoherence due to its interaction, represented by a quantum channel, with an environment. In this paper we study the dynamics of the mutual information $I(\rho_{AB})$, the classical correlations $C(\rho_{AB})$ and the quantum correlations $Q(\rho_{AB})$, as measured by the QD, in a two-qubit state the density matrix of which is the reduced density matrix obtained from the ground state of the TIM in 1d. The time evolution brought about by system-environment interactions is assumed to be Markovian in nature and the quantum channels considered are amplitude damping, bit-flip, phase-flip and bit-phase-flip. Each quantum channel is shown to be distinguished by a specific type of dynamics. In the case of the phase-flip channel, there is a finite time interval in which the quantum correlations are larger in magnitude than the classical correlations. For this channel as well as the bit-phase-flip channel, appropriate quantities associated with the dynamics of the correlations can be derived which signal the occurrence of a QPT.Comment: 8 pages, 7 figures, revtex4-1, version accepted for publication in Eur. Phys. J.